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Abstract 

This paper concerns the instantaneous frequency (IF) of continuous-time, zero-mean, complex- 
valued, proper, mean-square differentiable nonstationary Gaussian stochastic processes. We compute 
the probability density function for the IF for fixed time, which extends a result known for wide- 
sense stationary processes to nonstationary processes. For a fixed time the IF has either zero or infinite 
variance. For harmonizable processes we obtain as a byproduct that the mean of the IF, for fixed time, 
is the normalized first order frequency moment of the Wigner spectrum. 
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I. Introduction 



This paper treats the instantaneous frequency (IF) of mean square differentiable Gaussian zero-mean 
complex-valued proper nonstationary stochastic processes defined on R. The IF of a stochastic process 
is the derivative of the phase function and it is a stochastic process (real-valued) itself. Our main result 
is an explicit formula for the probability density function (pdf) for the IF for an arbitrary fixed time 
point. The pdf is parameterized by the covariance function of the process and its derivatives evaluated 
at the same time point. The time axis may be divided into three subsets with completely different IF 
process variance behavior: In the first subset, the pdf has heavy tails and behaves like x~ 3 for large x, 
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which means that the IF variance is infinite. In the second subset, the IF has a degenerate pdf consisting 
of a Dirac measure at its mean (that is, the IF is deterministic). In the third subset the IF is +00 with 
probability one. (Alternatively one may say that the IF is not defined on the third subset.) For wide-sense 
stationary (WSS) processes we show that the first subset is either R or empty, and in the latter case 
the covariance function has real part p x (t) = acos((3t) for a, (3 £ R, a > 0. For harmonizable but not 
WSS processes we show through examples that the first set may be R or it may be empty, but not ruling 
out other possibilities. The formula we obtain for the pdf of the IF for fixed time is a generalization of 
results on the pdf of the IF for Gaussian WSS processes derived by Miller [15] and Broman [2]. 

Our result is formulated for Gaussian zero-mean complex-valued proper processes that are mean- 
square differentiable and have mean-square continuous derivative. When we specialize to certain harmo- 
nizable stochastic processes z(i), the formula for the pdf of the IF for fixed time implies the identity 

where W z denotes the Wigner spectrum. The Wigner spectrum is the expected value of the Wigner 
distribution Wf, which is defined by 

W f (t, = / f(t + r/2)/(t - r/2)e- iT «dr, t,$£R, (2) 

for a function / : M h> C. 

Research about the IF has a long history in communications, signal processing and time-frequency 
analysis 0, H, ifTTl . In analog frequency modulation the IF (minus a constant carrier frequency) repre- 
sents the information in a modulated signal |fl7l . In time-frequency analysis [3] the IF of a deterministic 
signal / has been related to the Wigner distribution Wf by the formula 

* arg/w= raw (3) 

If a deterministic function is interpreted as a degenerate stochastic process, then this is a special case of 

CD. 

The Wigner distribution gives a time-frequency description of /. It was introduced in Quantum 
Mechanics as a candidate for a pdf of a particle in phase space. It is well known that Wf is only rarely a 
nonnegative function IfTTl so Wf may not be interpreted as an energy distribution or as a pdf in general. 
This is consistent with the Uncertainty Principle [ 1 1 1 which gives upper bounds on the resolution of phase 
space localization of particles, or, in the signal analysis interpretation, of the time-frequency resolution of 
signals. Nevertheless, if the Wigner distribution Wf is convolved by a sufficiently wide Gaussian function 
it becomes nonnegative, so domains of sufficiently large area in the time-frequency (phase) plane admits 
localization. The formula (01 supports the interpretation of Wf as a time-frequency distribution, since the 
right hand side is a normalized first-order frequency moment of the Wigner distribution, which delivers 
the center frequency of a narrowband signal. 
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The paper is organized as follows. After fixing some definitions and notation in Section [TT] we 
introduce the framework of mean-square differentiable Gaussian proper stochastic processes in Section 
Hill which also contains a background for the special case of harmonizable processes. In Section llVl we 
discuss and prove a precise version of the well-known formula (f3]) for the IF of a deterministic signal 
as a normalized first-order frequency moment of the Wigner distribution (for fixed time). In Section [V] 
we briefly give a background on earlier work on the IF of stochastic processes and its relation to the 
Wigner spectrum. Then, in Section |V|1 we prove our main result, which is a formula for the pdf of the 
IF for fixed time and its relation to the Wigner spectrum. Finally we show by examples in Section IVIII 
that a process may have infinite-variance IF for all time points or may have zero-variance IF for all time 
points. 

II. Preliminaries 

A probability space (CI, B, P) consists of a space CI, a cr-algebra B and a probability measure P. A 
random variable is a measurable function X : Cl H> R d (often denoted by a capital letter), and X induces 
a probability measure on B(R d ), the Borel cr-algebra, defined by P X {A) = Ppf" 1 ^)), A £ B(R d ). 
If the probability measure Px is absolutely continuous with respect to Lebesgue measure we have 
Px(A) = J A px(x)dx where px is the probability density function (pdf) of X. Sometimes, by abuse 
of notation, we will say that a probability measure that is a Dirac measure at a point a £ R d , denoted 
S a , has pdf S a . We denote the expectation of a random variable X by EX, and by Lq(CI) — Lq(CI, B, P) 
we understand the Hilbert space of second-order (finite variance) zero-mean complex-valued random 
variables. The space of second-order random variables with nonzero mean is denoted L 2 (Ct). 

Given a probability space (Cl, B, P), a continuous-time stochastic process on R is defined as a family 
of complex-valued S-measurable functions z t (lj) indexed by t £ R. We often suppress the variable u> £ CI 
and denote z(t) — z t (uj). Sometimes we write z t (oj) — z(t, lj) to emphasize the fact that z : R x Cl i->- C 
is a function of two variables. For fixed lj the function t H> z(t,Lj) is called a realization, trajectory or 
sample function. 

For a space U, we denote by \a the indicator function of the set Act/, that is xa(x) = 1 if 
x £ A and xa{x) = if x ^ A. The space of continuous functions on R is denoted C(R), the space 
of continuously differentiable functions on R is denoted (^(R), and the space of continuous functions 
decaying at infinity is denoted Co(R). This means that for any e > there is a compact set K £ C M. 
such that x E K £ =>■ \f(x)\ < e [21], where K £ = R \ K e denotes set complement. A derivative with 
respect to time is denoted / = df jdt. This notation is also used for stochastic processes where we define 
the derivative in the mean-square sense. That is, a process z(t) is mean-square differentiable at t = to 
if there exists z(to) £ L 2 (il) such that 

z(t + e) - z(t ) 



limE 



i(«o) 



0. (4) 
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For a function f(x, y) of two variables we write the partial derivative with respect to the first variable 
as dif(x,y) = df(x,y)/dx and with respect to the second variable as d2f(x,y) = df(x,y)/dy. The 
normalization of the Fourier transform for functions / G L 1 (R) used in this paper is 



&m = no = / me-^dx 

JWl 

which gives the inverse Fourier transform 

J JR 

We denote by t ^"L 1 (R) the space of functions / with Fourier transform / £ L 1 (M). For functions of 
several variables a partial Fourier transformation with respect to variable j is denoted ,^-j. The Wigner(- 
Ville) distribution (7), (9), ifTTI for / £ L 2 (R) is defined and denoted by (|2), where / denotes complex 
conjugate. 

A tempered distribution / belongs to the Sobolev space of order s£l (9), denoted / € H S (R), 
provided its Fourier transform / is locally integrable and satisfies 

(i + iei 2 n/(oi 2 ^<^. 



The Sobolev scale is a smoothness scale since for / € H S (M.) and s large, a certain amount of asymptotic 
decay at infinity of the Fourier transform is required. This implies that / will be differentiable to a degree 
that increases with s. 

Finally we use the notations M. + = [0, +00), R_ = (— 00, 0], the determinant of a square matrix M 
is det(M) = \M\ and the transpose of a vector x is x T . 

III. Mean square differentiable and harmonizable Gaussian stochastic 

PROCESSES 

Let z(t) = x(t)+iy(t) be a continuous-time, complex-valued, zero-mean Gaussian stochastic process 
z : R *-¥ Lq(Q.), not necessarily WSS. The assumption that z is Gaussian means the following J6). 
For any finite vector of time points t = (tj)™ =1 the sampled process real-valued 2n-vector Z t := 
(x(ti), ..., x(t n ), y(h), • • ■ , y(t n )) T has the pdf 

p Zt {u) = (27r)-' l |Afr 1 / 2 cxp l-]-u T Kr 1 ^\ , u € M 2n , 

provided the covariance matrix M = E,(Z t Zf) is invertible. A more general definition, which works 
also when M is singular, is the requirement that the sampled process Z t is a random variable with 
characteristic function 

<j>z t (u) = E(exp(m T Z t )) = exp ( u T Mu\ , u € R 2n . 
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The covariance function of the process z is denoted 

r z (t,s)=E(z{t)~z(sj), t,seR. 
We assume that z is proper lfl6l . that is E(z(t)z(s)) = 0, which implies 

r x (t,s)=r y (t,s), r yx (t,s) = -r yx (s,t), r z (t,s) =2r x (t,s) + 2ir yx (t,s), t,seR, (5) 

where r x (t, s) = E(x(t)x(s)), r y (t, s) = E(y(t)y(s)) and r yx (t, s) = E(y(t)x(s)). 

We require that the process z(t) be mean-square continuous and have a continuous mean-square 
derivative according to the following definition. 

Definition 1: A mean-square continuous process z(t) is mean-square differentiable with continuous 
derivative z(t) if there exists z(to) £ L 2 (Vt) such that (0|i is satisfied for all to £ R, and furthermore 
lim e _> E|i(t +e) - i(io)| 2 = for all t € M. 
This definition guarantees that d\r z , d\d2T z are continuous functions H, iflOll . lfl3ll and 

E(i(t)I(^)) = air 2 (t,s), E(z{t)~zjs))=d l d 2 r z {t,s), t,s£R. (6) 

Remark 2: Note that in Definition Q] we discuss a derivative process in the mean square sense only. 
That is, the definition does not imply that each realization is continuously differentiable with probability 
one. Conditions that are sufficient for continuously differentiable realizations with probability one are 
much more subtle and difficult iHJ, 0, lfT31 . 

For some results we will assume that z is strongly harmonizable (abbreviated in this paper as 
harmonizable) |[T2l . |[T3l . Il20l . which means that r z has a Fourier-Stieltjes representation 

r z (t,s)= (J e^- s ^m z (di,dri), (7) 



where m z is a measure of bounded variation on M 2 , called the spectral measure. This assumption implies 
that the process z(t) has a Fourier transform representation 

z{t) = I e^Z{d£), 



where the so-called spectral process Z : B(M) H> Lq(Q.) is a vector-valued measure of bounded 
semivariation |[T2l . 

An important special case of harmonizable processes are the mean-square continuous WSS processes. 
This means that there exists a continuous positive definite function p z such that r z (t, s) = p z (t— s) |fT3l . 
In the spectral domain the mean-square continuous WSS processes are characterized by m z (A, B) = 
H z (A n B), A, B £ B(R) for a non-negative bounded measure of one variable \i z . 

For harmonizable processes the following requirement in the spectral domain is sufficient to guarantee 
that the process z(t) is differentiable in the sense of Definition Q] 
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Definition 3: A harmonizable process z with spectral measure m z has spectral moments of order 
one if 

ff (1 + |£| 2 ) 1/2 (1 + \v\ 2 ) 1/2 \m z M,dr,) < co. (8) 

J JR 2 

Here \m z \ denotes the total variation measure of the complex measure m z 11211 . Note that this 
definition implies that all four integrals 

|ro,|(de,d»7), ff |^||»7||m,|(de,d»7), ff l£IK|(d£,cfy), ff \r)\\m z \(dt,dr,) 

J Jr 2 J Jr 2 J Jr 2 

are finite. Definition [3] guarantees that we may take partial derivatives under the integral in (O as 

d irz (t,s)= ff ite><X-"i)m,(dZ,dn), 



d 2 r z (t,s)= (-iv)e m -" Tl) m z (d^,dn), (9) 



did 2 r z (t,s)= // Zne^-^m^dri), 

J Jr 2 

due to Lebesgue's dominated convergence theorem ETj . Moreover, the process satisfies Definition Q] and 
© HH. 

The Wigner distribution (|2]i may be written Wf — &i{j ® / ° k) where j^jj denotes Fourier 
transformation in the second variable and k denotes the coordinate transformation 

K{x,y) = (x + y/2,x-y/2) «=>> K- l (x,y) = {{x + y)/2,x - y). (10) 

The Wigner (-Ville) spectrum [7| of a harmonizable process z(i) is defined by 

W, = &i(r t o k). (11) 

Since r z may not be an integrable function, the partial Fourier transformation in (fTTT i is in general defined 
with r z understood as a tempered distribution J9). However, in the case when r z 6 5o(K 2 ), which denotes 
Feichtinger's algebra ifTTl . we can write (fTTb as the partial Fourier integral 

W z (t,0 = f r z (t + T/2,t-T/2)e~ tTi dT, iel. (12) 

Under the same assumption (plus Gaussianity) we may interchange the order of integration and expec- 
tation so we also have 

W z (t, = E f f z(t + T/2)z(t - T/2)e- tTi dTJ =E(W»(t,0), *eM, 

that is, the Wigner spectrum is the expected value of the Wigner distribution of the stochastic process z 
defined by d2l El. 

Using (flOl l we may write the representation (0 as 

r z on{t,s)= ff e i ^ +t ^m z o K (d4,dri) (13) 

J Jr 2 

and thus by identification with (fTTT > it follows that 

W z (t, d£) = 2tt f e ltv m z o *(#, drj). (14) 
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IV. The IF and the Wigner distribution 

The argument (or phase) of a complex number z = x + iy eC\ {0} is defined by z = |z|e largz 
where we impose the restriction — 7r < arg z < ir. We have 

arctan (^) , x > 0, 

7rsgn(y) + arctan (f) , x < 0, y ^ 0, 
air, = < |sgn(y), x = 0, y^O, (15) 

7T, X < 0, y = 0, 

undefined, .t = y = 0. 

For z G U := C\{(— 00, 0] +i0} we have arg z = Im (logz) = —i\og(z/\z\) where the principal branch 
of the logarithm function is understood, i.e. — 7r < Im (log z) < n. The function logz is holomorphic 
from the domain U onto R + i(— 7r,7r) C C, and therefore z H> argz is a smooth function U i-> (— tt,tv). 
Let fit) = x(t) + iy(t) be a function / : R 1-4 C If /(£) ^Owe reserve the notation 

¥>(*)= arg /(t) 

for the phase function. For / G C(R) the set L// := {t G R : f(t) G [/} C R is open and ^ = argo/ : 
[// n- (— 7T, 7r) is continuous . If, moreover, / G C 1 (R) then it follows from above that ip is differentiable 
on Uf, and since 4? arctan(i) = 1/(1 + t 2 ) the derivative is 

■ M _ d ^ > _ 1 x(t)y(t) - ±{t)y(t) _ x(t)y(t) - x(t)y(t) 

at I + tttt x ^> x(t)+y(t) 

for t G Uf. In fact, for {(eM: x(t) > 0} and {t : x(t) < 0, y{t) ^ 0} this follows from (Q3)- In the 
remaining case {t € R : x(£) = 0, y(t) 7^ 0} we may use the following modified definition, equivalent 
to (IB) in {t G R : y(t) ^ 0}, 

7T , . ( X 

argz = — sgn(y) — arctan — 

2 \y 

which gives (pit) = —x(t)/y(t) for {t G R : x(£) = 0, y(t) 7^ 0}. This also leads to the expression on 
the right hand side of ( TT6J I. 

Clearly the formula ( PToT ) can be extended from the domain Uf to {t : /(i) 7^ 0} and p is still 
continuous {t : /(i) 7^ 0} h» R. The instantaneous frequency (IF) (5), J4|, G) of /(£) is defined by 
( [ToT l as the derivative <£(£) with domain {t : f(t) 7^ 0}. For an exponential function e^' with frequency 
£ the IF is thus £ constantly, which means that the term instantaneous frequency is an extension of the 
concept of a constant (global) frequency. For {t : f(t) = 0} it will turn out to be convenient to define 
(by abuse of notation) ip(t) = +00 so in summary we have for / € C 1 (R) 

{ x{t)y(t)-x(t)y(t) , f 2 (f\ , ..2( f \ ^ n 

*w&m lf x w + vw>°. (17) 

+00 if x 2 (t)+y 2 (t) = 0. 

If / is real-valued and continuous then f(t) 7^ implies that / has constant sign, that is ip(t) = or 
<p(t) = tt, in a neighborhood of t, and hence <£>(£) is well-defined and equals zero in this neighborhood 
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of t. The derivative of the phase function of real-valued signals is thus not very interesting. However, the 
Hilbert transform J4|, ifTTl of a real-valued signal gives rise to a complex-valued signal with a nonzero 
IF. Since this transformation from a real- valued signal to a so-called analytic signal transforms cos(£i) 
into e'**, it gives a natural definition of the IF of a real-valued signal, commonly used in the literature. 

There is a connection between the IF of a sufficiently smooth and decaying function and the Wigner 
distribution 0, Q. The heuristic version of this result is well known [3] but here we prove a more 
precise statement. 

Proposition 4: Suppose e > 0, / G Hi +e (R), f(t) = x(t) + iy(t) and cp(t) is defined by (ITTb - 
Then for any tgl such that f(t) ^ we have 

J R W> (*>£)<£ 

Proof: The assumption is 

/(i + iei 2 ) l+e i/(0l 2 ^<oo. d9) 

Jm. 



Let t G R be fixed and arbitrary such that /(*) ^ 0. Define the function g(r) := f(t + r/2)f(t - r/2). 
The requirement (fT9l implies that / G L 2 (R) which means that / G L 2 (R) by Plancherel's theorem. 
Thus by the Cauchy-Schwarz inequality we have 



g = f(t + -/2)f(t--/2)€L l (R). (20) 

The Cauchy-Schwarz inequality and ( fT9] l give 

(i + iei 2 ) 1/2 i/(e)R = Jq. + i£i 2 ) |+i+| -*- f l/(« 

1/9 1/9 \^ *- ) 

< Qf a + iei 2 ) l+E i/(oi 2 ^) (jT(i + iei 2 )-*- s ^) < oo. 

Hence / G L^M) and since &(f{t + -/2))(£) = 2e l2 ^/(2£) we have J?(/(i + -/2)) G L 1 ^)- Thus 



&{f(t - -/2)) = <£"(/(< + -/2)) implies that 



^.g = &(f(t + -/2)/(t - -/2)) = — ^(/(t + -/2)) * ^(/(* + 72)) G L\R). (22) 

Moreover, / G L 1 (R) implies that / G Co(R) according to the Riemann-Lebesgue lemma. Also, since 
J^"/(£) = «£/(£) (where / is considered a tempered distribution), and £/(£) G L 1 (R) according to 
OTT i. we have / G Co(R) again by the Riemann-Lebesgue lemma. Thus by (|20| > and (l22l we have 
g G L 1 (R) n t ^"L 1 (R) n Co(R), which means that Fourier's inversion formula holds pointwise |21| for 
the function g: We have g(r) = (jF -1 g)(r) for all r G R. In particular 

2n\f(t)\ 2 = 2^(0) = f g($)d£ = / W f &£)<%, (23) 

Jw Jw 

where the last equality uses the definition (|2}. 
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Concerning the numerator in (fl8b we use integration by parts and the just proven fact that / vanishes 
at infinity to obtain 

d 

{e~°"»)dT 

(24) 



eW>(t,0= / f{t + T/2)f{t-T/2)ij-{e-^) dT 



d_ 



< ' --(/(t + r/2)/(t-r/2))dr. 



Let us study the function gi(r) = f(t + r/2)/(i — r/2) and 

d 



h (r) = Tt (/(< + r/2)/(t - r/2)) = ±( 5l (r) - 5l (-r)). 

The assumption (O implies that £/(£) € L 2 (R) and thus J*"/ € L 2 (R), and / e L 2 (R) by Plancherel's 
theorem. Since we already know that / e L 2 (M.) the Cauchy-Schwarz inequality gives g\ £ i 1 (R), 



From above we know that &(f(t- -/2)) = J*"(/(i + -/2)) = 2e _ * 2£ V(2£) € I/ X (R). Likewise we 
have J?(/(i + -/2))(£) = 2e l2 «*/(2£) = 2ie l2 «*2£/(2£). This function belongs to ^(R) because of 
(EB- Thus 

?i = ^(/(* + -/2)) * JW--/2)) G ^(R). 

Hence we have proved that g^ e Co(R)nL 1 (R)n^L 1 (R) which means that Fourier's inversion formula 
holds for gi and hence also for h. Denoting f(t) = x(t) + iy(t), integration of (l24l gives 

ZW f {t,OdZ = -%( &h{Od£ = -2mh(0) = -m (f(t)W- /(*)/(?) 

Jr v / (25) 

= 27r(x(t)y{t) - *(*)»(*))• 
Finally < fT~8l > follows from a combination of ( fT71 ), d23t and d25l i. ■ 

V. The IF for stochastic processes and the Wigner spectrum 

One of our goals in this paper is to generalize Proposition [4] from deterministic functions to certain 
stochastic processes defined on R, denoted z{t) = x{t) + iy(t) = \z{t)\e lv ^ . More precisely we would 
like to prove the formula 

E(0(t)) = J J^LJ12 fi , (26) 

In fact, we will compute the pdf of the random variable <fi(t), for fixed t € R, and then as a consequence 
derive the formula above. This problem has been studied by Miller ifTSl and Broman ||2) for WSS 
Gaussian proper stochastic processes. Miller and Broman independently derived the probability density 
function for (p(t) for fixed t, using either the hypothesis that the process is proper (in lfl5l ) or the more 
restrictive hypothesis that the signal is analytic ifTTl (in 0). (See also lfl9l .) Our aim is to generalize 
their results from WSS to Gaussian proper nonstationary processes that have mean-square continuous 
derivative. As a special case we will study certain harmonizable processes and prove the identity ( |261 ) 
for them. 
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For processes more general than WSS it has been customary, in parts of the literature, to define the 
instantaneous frequency as Q, lfT4l 

0(();=£ (iM|g*M, mtr>0 . (27) 

This definition replaces the random variable in the denominator of dTvT > by its expected value. It gives 
a generalization of the formula ( fT~8T > of Proposition [4] with the IF <p(t) replaced by E(ip(t)) and Wf 
replaced by Wf, as follows. 

Proposition 5: Suppose z(t) is a proper harmonizable process such that the spectral measure satisfies 
©. Then for any t e R such that E|z(i)| 2 > we have 

Ett(t)) = Hm ^ — -. (28) 

Proof: First we note that (0, ( TTOb and ( TBI give 

E\z(t)\ 2 = r z (t,t)= ff e u ^-^m z {di,dr 1 )= ff e^m z o «(d£, drf) = ±- f W z {t,d£). (29) 

J JR2 JJR2 /7T J^gR 

Next, since E(x(t)y(t)) = dir yx (t,t) and E(x(t)y(t)) = d 2 r yx (t,t), we obtain using ©, © and (O 
E (x(i)y(t) - i(t)|/(t)) = dir vx (t, t) - d 2 r yx {t, t) 

= ^(d x r z (t,t)-d 2 r z (t,t)) 

=H[b { s +v) ^ )mM > dr)) ) (30) 



£,e Un m z on(d£„dri) 

yJ ZW z (t,dO. 

Z7r JfeR 



Finally <J28]» follows from a combination of d27]), d29j) and (O. ■ 

The definition d27b is however not the natural definition of the IF stochastic process. Instead we 
define the IF by < fTvT > where z(t) = x(t) + iy(t) = \z(t)\e lLpi ^ is a zero-mean, complex-valued, proper, 
Gaussian process that is differentiable in the sense of Definition Q] Thus ip(t) is a IU {+oo}-valued 
stochastic process defined on M. 

Remark 6: A trivial case of a second-order zero-mean stochastic process consists of a random 
variable X £ L 2 (Vl) times a function / e C 1 (R), that is z(t) — Xf(t). For such a process we have 
z(t,uj) = \X(u)\\f(t)\e i ^ TsX ^ +v ^\ where f(t) = \f(t)\e lv(t \ provided /(*) ^ and X(u) ^ 0. It 
follows that the IF of z{t,u) is ip(t) provided that f(t) ^ and X(u) ^ 0. That is, the IF of z(t) is 
t P{t)X{u:X(ui)yio}{ UJ ) + °°X{uj:X (uj)=o} {u) ■ This means that the IF is essentially the deterministic IF of / 
and is only stochastic in the sense that for certain u € CI it is +oo for all igf, and for the remaining 
lo it does not depend on ui. This is true also when X is non-Gaussian and improper, that is EX 2 ^ 0. 
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VI. The probability density function for the stochastic IF FOR fixed time 

In this section we will derive the pdf of the IF stochastic process ip(t), defined by (TPTb for a mean- 
square differentiable stochastic process z(t) = x(t) + iy(t) = \z(t)\e l(p ^\ for a fixed arbitrary t G E. As 
a consequence we will obtain the formula d2o*l ) for E(ip(t)) provided the process is harmonizable. 

Let z(t) = x(t) + iy(t) be a proper Gaussian stochastic process which is differentiable according to 
Definition [T] Fix t£l and define the random column four-vector X = (x(t), y(t),y(t),x(t)) T . Then X 
is zero-mean and Gaussian, since Gaussianity is preserved under mean square limits |6|. Because z(t) 
is proper, © holds, and thus Ex(t)y(t) = Ey(t)x(t) = d x d 2 r yx (t, t) = and Ey(t)y(t) = dir y (t, t) = 
dir x (t,t) = Ex(t)x(t). Hence X has covariance matrix M — EXX T G R 4x4 with the structure 



/ 



M = 



a 


6 





c 


b 


d 


c 








c 


a 


-6 


c 





-b 


d 



\ 



(31) 



and parameter values 



a = Ex(t) 2 = r x (t,t), b = Ex(t)y(t) = dir vx (t,t), 

c = Ey(t)y(t) = dir x (t,t), d = Ey(t)y(t) = d 1 d 2 r x (t,t). 



We have \M\ = (ad- 



b 2 ) 2 and 



ad 



b 2 >0, 



(32) 



(33) 



which follows from taking the determinant of the upper left 3x3 submatrix of the nonnegative definite 
matrix M. 

Our main technical result concerns real Gaussian zero-mean four-vectors with covariance matrix < T3TT >- 
In the proof we will need a small lemma that may be considered "folkloristic" in probability theory and 
not necessary to prove. Nevertheless we prove it for completeness and clarity. 

Lemma 7: Let X : Q h^ W l be a random variable. Suppose that N e B(R d ), x e R d are two 
disjoint null sets for P x , i.e. P X (N) = P x ({x }) = 0, x £ N. If we define X : Q h-> R d \ N by 

X(u), LuenXX-^N), 

' \ v j, (34) 

^ xo, weX-^N), 
then P X (A) = P X (A) for all A e B(R d ). 

Proof: Let A e B(R d ). We write A = A x U A 2 U A 3 as a pairwise disjoint union with Ai = A n iV, 
A 2 = (An {a; }) \ A^ and A 3 = A \ (N U {xo}). Since X takes values in M d \ A^ we have P X (A X ) < 
P X (N) = 0, and moreover P X (A 2 ) < P x ({x a }) = P(X- 1 (N)UX- 1 (x )) < P X (N)+P x ({x }) = 0. 
Thus P X (A) - P X (A X ) + P X (A 2 ) + P X (A 3 ) =P X (A\(NU {x })) = P X (A\(NU {x })), where 
the final equality follows from X(u) G A \ (N U {x a }) <=> w G X _1 (^) \ (^ _1 (^) U X-^Zo)) <=> 



x( w ) = 
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X(u) EA\{NU {xo}). This finally gives P X (A) = P x (An(NU {x })) + P X {A\(NU {x })) = 

P x (A\(NU{x })) = P X (A) because P X (A n (N U {x })) < P X (N) + P X ({x }) =0. ■ 

Next follows the crucial result about Gaussian four-vectors with covariance matrix M according to 

(ED- 

Proposition 8: Let X — (X\, X2, X3, X4) be a zero-mean Gaussian real-valued vector with covari- 
ance matrix M defined in OTI ). Define the R U {+oo}-valued random variable 

Xl x)+xf* ' if X l + X l > °> 
+00 if xf + xf = 0. 



F 



(35) 



If ad — c 2 — b 2 > then Y has probability density function 

Pr{y) = l\M\ 1 ' 2 ((ay -b? + \M\W 
Consequently Y has infinite variance and mean 

EY 



-3/2 



b 
a 



(36) 



(37) 



If ad 



we have the two subcases: 



M- 



\M\ 



-1/2 



-c\ 



V 



(38) 



(i) If a > then Y = b/a a.s. 
(ii) If a = then Y = +00 a.s. 
Proof: Suppose first that ad — c 2 — b 2 > which means that M is invertible with inverse 

f d -b 

—b a —c 

-c d b 

—c 6 a / 

Furthermore ad—c 2 — b 2 > implies a > 0. Define the smooth function / : M + x (— 71", 7r] xKxl 1 — ► M 4 

h(y) = yi cos 2/2, 

/a(y) = y3siny 2 +2/12/4 cos y 2 

/s(j/) = 2/i sin 2/2, 

/4O/) = 2/3 cos 2/2 - 2/12/4 sin 2/2 



(39) 



The Jacobian of / is 



Df(y) 



( 



cos 2/2 -2/1 sin 2/2 

2/4 cos 2/2 2/3 cos 2/2-2/12/4 sin 2/2 sin 2/2 2/1 cos j/ 2 
sin 2/2 2/1 cos 2/2 



\ 



V - 2/4 sin 2/2 - 2/3 sin 2/2 - 2/1 2/4 cos 2/2 cos 2/2 - 2/i sin 2/2 / 



July 8, 2010 



DRAFT 



13 



whose determinant is det Df(y) = y\. Note that (2/1,2/2) H> (fi{y), h{y)) is the polar-to-rectangular 
coordinate transformation on R 2 . 

We will use / as a coordinate transformation, and then we will need / to be a bijection with a 
differentiable inverse. Since /(0, 2/2,0, 2/4) = for any 2/2 € (— 7r, 7r] and any 2/4 £ R, the function / 
is not injective on the domain R + x (— tt,tt] xlxl. Therefore we need to restrict the domain of /. 
Define 

N' = ({0} x (-7T, irjxlxl) (J(R+ x {tt} x R x R) c R+ x (-tt, tt] x R x R (40) 

and 

JV = {iel 4 :ii<0,i3 = 0}=i-xMx{0}xlcl 4 . (41) 

Then it can be verified that f(N') — N and the restriction 

/ : R + x (-tt, vr] x R x R \ iV' 1 — > R 4 \ N (42) 

of / to the open set R + x (— n, tt] x R x R \ N' is surjective and injective. Its inverse is / _1 = g = 
(91, 92, 93, gi) where 

fll(s) = (*? + *§) 1/2 , 

<7 2 (a;) = arg(xi +ix 3 ), 

X2X3 + X1X4 (43) 



93 (x) 



(xi+xiy/*-' 

X\X 2 - X3X4 



9i(x) - 

Since x ^ ./V => a; 2 + a; 2 > it is clear that gx, 2/3, and 2/4 are differentiable on R 4 \ N. Since the 
nonpositive xi-axis in the (xx, £3)-plane, in conjunction with any (#2) £4) € R 2 , does not belong to 
R 4 \ N, 272 is differentiable on R 4 \ N. Hence the restriction of / (l42l is differentiable and has a 
differentiable inverse. 

The sets N' C R+ x (—tt, tt] x R x R and iVcM 4 are null sets with respect to Lebesgue measure, 
and Px{N) = 0. Let x e R 4 \ N be fixed arbitrary and define the random variable X : tt !->■ R 4 \ N 
by 

%) j w «»u-w (44) 

[ zo w€X- 1 (AT). 

It follows from Lemma|7]that X and X induce identical probability measures, i.e. P X (A) — Px(A) for 
all A e 2?(R 4 ), which means that 

Pjf (j4) = (27r)" 2 |Afr 1/2 / exp (-^x T M~ l x\ dx, A e £(R 4 ), (45) 

because X is a zero-mean Gaussian with covariance matrix M, 
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Define the random variable Z : O H> R 4 by Z = / _1 (X), which is well defined because X(fi) = 
R \X and the restriction d42l i is bijective. If we write Z = [Z\, Z%, Z3, Z4) and xq = (xo,i, £0,2, ^0.3, ^0.4) 
then according to ([35), d43l and d44l) we have 

_, , f ^M. uen\x-\N), 

Z a {uj) = { (46) 

— 3.2 +a ,2 , w e a (jv ;. 

V ^"0,1^ 0,3 

We claim that the random variables Y and Z4 have identical probability measures, that is 

P Y {A) = P Zi {A), AeB{R). (47) 

In fact, let A G B(R). We decompose Y- l {A) = {Y- l {A) n X _1 (X)) UO^O 1 ) \^ _1 (A0). Since 
P^" 1 ^) n X-^AT)) < PfX-^X)) = P X (N) = and Y- l {A) \ X-^N) = Z^(A) \ X~ 1 (A r ) 
according to(|46), we obtain iV (A) =P(F- 1 (^)) = P(y- 1 (A)\X- 1 (A^)) = P(Z^ 1 (,4)\X- 1 (X)) = 
P(Z^V)) = P Zi (A)> proving g3. 

Denote the probability density functions for X, X and Z by p^, p^ and pz, respectively. Then 
d45l > says that p%{x) = Px(%) = (27r)~ 2 |A/| _1/2 exp (-^x T M~ 1 x), x G K 4 . For an arbitrary Borel 
set A G B(R + x (-7T, tt] x R x R) we have 

Pz(A) = P(Z G A) = F(Z G i4 \ (X G X)) = P(X G /(A) \ (X G X)) 

= P(X G /(A) \ X) = P(X G /(A \ X')) 



p x (x)dx = p x o f(y)\detDf(y)\dy 

f(A\N>) JA\N> 

Pxo f(y)\detDf(y)\dy 



= / Px ° f(y)\detDf(y)\ X[o,+oc)(yi)X(-Tr,Tr](y2) dy. 

J A 

In fact, the seventh equality above is the formula for changing variables in integrals [8, Theorem 5.8]. 
To justify its use, we need the fact that (1421 is differentiable and has a differentiable inverse, which has 
been proved above. On inserting d38l and (139b we obtain 

Pz{z) =Px(/(2 ; ))|det£)/(z)|x[o,+oo)(2i)X(_^, 7r ](22) 

= (27r)- 2 |M|-^ 1 2 exp \-\f{z) T M- x f[z)\ X[o,+oo)(*i)X(-^](^) 

= (27r)- 2 |A/r^z 2 exp (— |Af|~4 (z 2 (d + az| - 26z 4 ) + azf - 2cziz 3 ) /2) X[0,+oo)(«i)X(-7r, w ](^) 

(48) 

after some computations. Using J °° x 2 exp(— sx 2 /2)dx = y / n/2s~ 3 / 2 , s > 0, and |M| = (ad—c 2 — b 2 ) 2 , 
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we obtain the marginal probability density for Z± 

/>00 />7T />00 

Pz i { z i)= III Pz(zi,Z2,Z3,z 4 )dz 1 dz 2 dz 3 



= (2tt)- 1 \M\- 1 ' 2 / z 1 2 exp(-|Af|- 1 / 2 (z 1 2 (d + az2-26z 4 )/2)) 

x ( / cxp (-\M T 1/2 (azl - 2cziz 3 ) /2\ dz 3 \ dz x 
= (27ra)- 1 / 2 |Af|- 1 / 4 / z 2 exp (-|M r 1/2 z 2 (d + az 2 - 26z 4 - c 2 /a)/2) dzi 

= a -\M\^([a Zi -bf + \M\^y V2 - 
The earlier noted observation J47l ) now proves d36l ). Finally, ^(x(x 2 + s)" 1 / 2 ) = s(x 2 + s)~ 3 / 2 gives 



,■2 



(a;(x 2 + s) _1/2 )dai = -, s > 0, 



l.r- + s) 3 / 2 S 7 ^ 



so we have 



EZ 4 = ||M| 1 /2 if 



zdz 



oo ((o^-^a + IMIVa) 3 / 2 



2a ./-oo (z 2 + |M| 1 /2) 3 / 2 
6 



a 

,-3 



Note that pz A (y) behaves like Cy for large y. Therefore 

(y - EY4) 2 pzAy) d y = +°° 



which means that the variance of Z4, and therefore also that of Y, is infinite. This proves the proposition 
in the case ad— c 2 — b 2 > 0. 

It remains to consider the case when ad—c 2 — b 2 — 0, i.e. M is not invertible. If a = then Xi = 
a.s. and X3 = a.s. =>■ y = +00 a.s. This proves case (ii). Assume henceforth that a > 0. If a > and 
<2 = then X 2 = a.s., X 4 = a.s. and X 2 + X 2 > a.s. =i- y = a.s. Since 6 2 + c 2 = arf = we 
have & = so case (i) is proved if d = 0. 

The rest of the proof is devoted to the case (i) with ad—c 2 — b 2 = 0, a > and d > 0. We will employ 
a regularization technique. The characteristic polynomial of M is det(A7 — M) = A 2 (A — (a + d)) 2 . Let 

Ua+d = M{M -( a + d)I)C M 4 

denote the two-dimensional eigenspace for M corresponding to the nonzero eigenvalue A = a + d. The 
probability measure Px is a Gaussian which is supported on U a +d and non-degenerate on this two- 
dimensional subspace. With N defined by ( |4TT >, it can be verified that N n U a +d C W C R 4 where VF 
is a linear subspace with dimW^ = 1. This means that Px(N) = 0. If X is defined by (|44| | for some 
£0 e M 4 \ iV, then Lemma [7] gives P x = Pg. 
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Let n > be an integer and define iW = (X 1 ,X 2 + X' 2 /^/n, X 3 , X 4 + X'Jy/n) where X' 2 , X' A 
are Gaussian zero-mean, unit- variance random variables, pairwise independent of each other and of X\, 
X 2 , X3 and X4. Then X 1 - 11 ^ : h-> M 4 is a Gaussian random variable with covariance matrix 

1 a b c \ 

6 d+i c 

Oca —6 

c -b d+i 

the determinant of which is |M„| = a 2 /n 2 > 0, so M„ is invertible. For the characteristic functions of 
X^ 71 ' and X, denoted by <f>xto and cf>x respectively, we have 



M„ 



V 



/ 



x{ „)(O=E(exp(^ T xW)) 

1 



CX P [ -nfMnt 



exp 



■^e-£(d+& 



exp(-~^M£) = <M£)> 



— > +00, £ e w. 

This is equivalent to the weak convergence of probability measures PxM —* Px as n 



-00 



II 



l( n »(w) 



Theorem 7.6], which means that J h(x)P X (.n) (dx) — > J h(x)Px(dx) for all bounded and continuous 
functions h defined on M 4 . 

If we define the random variable X ( n ) by 

XW(«), we!l\(XW)- 1 (iV), 

x w^lW)" 1 ^). 

with A^ defined by (l4TT i and £0 G M 4 \ A" fixed for all n > 0, then Lemma [7] again gives Py(»o = ^x(«) 
for all n > 0, so we have P^ (n) -> Pj weakly. Let us define Z = / _1 (X) and #(") = /-'(fW) 
where Z = (Z 1 ,Z 2 ,Z 3 ,Z A ) and #W = (z[ n) , Z { 2 \ zf> , z[ n) ). Then we have P Y = P? 4 as in the 
first part of the proof. Since / _1 is continuous on the range spaces of X and X^ n \ and Pj(„) 
weakly, we may conclude that P^w — > Pz weakly |fl] Section 1.5]. The weak convergence P^m 
is equivalent to the limit of characteristic functions [lj 

n— >oo 

In the following we will compute <f) Z (.n) (£) using the probability density computed in (l48l l for invertible 
covariance matrix M. Since | TA^^I 4 / 2 = a/n we obtain from (l48b the pdf of #(") as 

p Z (n)(«) = a(27] _, 2 zicxp(^-— (zi(d+ l/n + azj -2bz i ) + azl - 2cz 1 z 3 )j X[o,+oo)(zi)X(-*,n](z2)- 

(50) 
It depends trivially on z 2 , so we may concentrate on the M 3 -valued random variable U^ n > = (Z{ , Z\ , Z4 ). 
The marginal probability density for LA") is 



Pz 



(49) 



Py(n)(2l,«3,«4) = / p Z (.n)(Zl,Z2,Z3,Z4,)dZ2 

zlexp\-— (z 2 (d + l/n + azl- 2bz 4 ) + azl - 2cZiZ 3 )J X[o,+oo)(^i)- 



— 7T 

n 



2ira 
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Thus the characteristic function of U^ n ' is 

&/<»)(£) =<hw (&,&,&) =E(exp(i(^ ) 6,^) T C/ ( " ) )) 

Pt/(n) (zi, z 3 , z 4 ) exp(?(zi^i + z 3 £ 3 + z 4 £ 4 ))<izi<iz3<i2:4 

+°° '•+ 00 / n o / 2 2b \ n ( 2 2c 

: 4 ~ ^4 ) + ^ 4 £ 4 - — ( Z 3 Z ZlZ3 I ■ ' " ;C; I " ■ '*''' ' i J " - J 

— 00 J — 00 



exp ( --^1 f z\ z 4 j + iz 4 £ 4 - - I z\ ~-Z\Z$ ) + iz 3 £a ) dz 3 dz 4 ) dz\ 



( b &\ f ( z ( c \ £ \ 

= a" 1 exp U£ 4 ---L / 2l exp --i + izi ^i+^ 3 - -— ^ d^i, (fi,&,£4) e 

\ a 2ny J \ 2a V a/ znzfj 

Set [/ = (Zi, Z 3 , Z 4 ). It follows from ( 149} and dominated convergence that 

<M£i,£3)60 = 0z(£i> o )£3>&) = lim <Asc»)(£i» 0,6,60 = Hm 0cj(»)(Ci.6)^4) 

n ^00 n »-oo 

exp(i£ 4 -W ziexp f --i- +izi ^1 + £3- J J d^i, (€i,£3,60 G 



a- 1 



(51) 
From (T5Tb we conclude that the characteristic function for Z 4 is 

4>zAU) = 0f/(0, 0,£ 4 ) = a -1 exp ( ?£ 4 - j / ziexp f --i- j dz\ = exp f i£ 4 - j , 

which implies that the probability measure for Z 4 is P% A = 8 b / a . Finally this gives Py = Pzi = <V a =>■ 
Y = b/a a.s. ■ 

Remark 9: If z is WSS then r x (t, s) = p x (t—s) for some even function p x , which implies dp x (0) = 
=>■ c — dir x (t, t) = 0. In this case the matrix M is block-diagonal and this special case of Proposition 
[8] was proved by Miller [15]. 

As a corollary to Proposition [S] we obtain the following result for Gaussian stochastic processes. To 
formulate it we need to introduce a partition of the time axis into two disjoint sets, depending on the 
covariance function r z for a given process z. The partition is 

T:={teR: r x {t,t)d x d?,r x {t,t) - {d iryx {t,t)) 2 - (dir*(M)) 2 = 0}, 

(52) 
=>T c = R\T = {teR: r,(M)<9icV x (M) - (dir yx (t,t)) 2 - (d x r x {t,t)) 2 > 0}, 

since r x (t, t)d x d 2 r x (t, t) - (d x r yx (t, t)) 2 - (d x r x (t, t)) 2 > holds for all t e K due to H3 and 0D- It 
is clear that TCRis closed since it is the inverse image of {0} of a continuous function. We further 
subdivide T = T" U T", where T" C R is closed, as a disjoint union of the measureable sets defined by 

T' :={teT:r x {t 7 t)>0}, 

T" :={teT:r x (t,t) = 0}. 

Corollary 10: Let zit) = xit) + iy(t) = \z{t)\e lip ^ be a zero-mean proper Gaussian stochastic 
process which is differentiable according to Definition Q] and let the instantaneous frequency stochastic 
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process ip(t) be defined by (fTTI i for all tcR. Fix tel, define a, b, c, d by (f32l > and the matrix M by 
(f3Tb . Then tp(t) has pdf 

/ ||7Wf|V2( (ay _ 6) 2 + [M |l/ 2 )-3/2 [fteR \ T: 

Pv(t)(y)=< (53) 

I h/ a (y) if t€T', 

and if i E T" then y>(t) = +00 a.s. Consequently (p(t) has mean 

[ +00 if i e T", 

and variance 

+00 if t € R \ T, 

:i;-(/ 1 -e:,-(/d- = <; o ifter, 

undefined if (ef. 

Finally we restrict to harmonizable processes whose spectral measure satisfies ||8), which admits a 
connection to the Wigner spectrum as follows. 

Theorem 11: Let z{t) = x(t) + iy(t) = \z{t)\e lv ^ be a proper Gaussian harmonizable stochastic 
process whose spectral measure satisfies (H) and let the instantaneous frequency stochastic process ip(t) 
be defined by O for all t € M. Fix i e M. Then <£(*) has pdf d53| where a, 6, c, d is defined by 02) 
and the matrix M by PTT l. The process (p(t) has mean 



I ef *ew,(t,dt) 



if £ e R \ T", 



E0(*) = ^ J*« W *M« 

-00 if £ e T", 

and variance 

-00 if t € R \ T, 

EM./!- Ei(/.).)' J = { ifteT', 

undefined if £ e T". 

Proof: First we observe that (0, (|29T l and d32l give 

a = r x (t,t) = - Re r*(i,i) = -r,(t,t) = — / W,{t,dg). 

Next, from the requirement I© it follows that (fl~3T > may be differentiated with respect to r under the 

integral: 

£ (r, o «) (i, r) = // z£e^« +t ">m 2 o *(#, A;). 
or JJ R 2 

Using (TBI , this gives 

£ (r 2 o «) (i, 0) = i ff &»>m, o «(#, drf) = ^- ( £W z (t, d£). (54) 

OT JJ R 2 27T J ?eR 
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On the other hand (0 implies d\r x {t,t) — d 2 r x (t,t) and d2r yx (t,t) = —dir yx (t,t), which gives 



d 

— (r z o k) (t, t) 



d 

= 2— (r x o k + ir yx o k) (t, t) 

t=0 "T 



T=0 

(55) 



= dir x (t,t) -d 2 r x (t,t) +i(d 1 r yx (t,t) - d 2 r yx {t,t)) 

= 2id\r yx {t,t). 
Combining d54l >, ( l55l > and ( |32l we have proved 



4tt& = f £W z (t,dO. 

JfGR 



The result now follows from Corollary [10] which in particular gives 



a I^W z (t,dO 



provided tel\f. 



VII. Cases of constantly zero or constantly infinite variance IF 

The process ip(t) exhibits completely different behavior on t <G T" = T\T" (where it is equal to b/a 
with probability one), and t £ M. \ T (where it has infinite variance with mean value b/a). It is therefore 
of interest to investigate questions like necessary or sufficient conditions for T = or T = R. We restrict 
to harmonizable processes whose spectral measure satisfies (0. First we look at WSS processes. 

Proposition 12: Suppose z satisfies the requirements of Theorem fTTI and z is WSS with covariance 
function r z (t, s) = p z (t — s) and nonzero. Then we have T" — 0, and either T = or T = T" = R. In 
the latter case p z — 2p x + 2ip yx where p x (t) = acos((3t), a > and /? > 0. 

Proof: The assumption that z is nonzero means that p x (0) > 0. As observed in Remark [9] we have 
c = dir x (t, t) = if z is WSS. Since d\&2r x (t, s) = —d 2 p x (t — s) we have 

r = {t£l: r x {t,t)d x d 2 r x {t,t) - (dir>(£, t)) 2 } = {t £ R : -Px(0)<9 2 ^(0) = (d/^O)) 2 }. 

It is thus clear that either T = R or T = 0, depending on whether — p x (0)d 2 p x (0) — (dp yx (0)) 2 or 
-p x (0)d 2 p x (0) > (dp yx (0)) 2 is satisfied. If T = then of course T" = 0. Thus suppose that T = R, 
i.e. — p :E (0)9 2 p 2 ;(0) = (9yO ya; (0)) 2 . This means that 

Ea;(t)x(t)Ex(t)x(t) = (Ey(t)x(i)) 2 < Ex(t)x(t)Ey(t)y(t) = Ex(t)x(t)Ex(t)x(t), t £ R, 

where the inequality is the Cauchy-Schwarz inequality, and where Ex(t)x(t) = Ey(t)y(t) = d\d 2 r x (t, t) 
since r x = r y . Therefore we have equality in Cauchy-Schwarz, Ex(t)x(t)Ey(t)y(t) = (Ey(t)x(t)) 2 for 
all t€l. This means that y(t) equals a real multiple of x (t) as a member of Lq(O), that is y(t) = c(t)x(t) 
where c(t) £ R. This gives r y (t,s) = did2r y (t,s) = —d 2 p x (t — s) = c(t)c(s)p x (t — s) => c(t) = 



±y— 9 2 / o a: (0)/p a ;(0), that is c(i) is constant. We obtain the differential equation —d 2 p x (t) = c 2 p x (t) 



July 8, 2010 DRAFT 



20 



with solution p x (t) — acos(ct) = acos(\c\t) (since p x is even) where a = p x (0) = Ex(t) 2 > 0, and 
P = \c\ > 0. Finally, T = T', that is T" = 0, since ^(0) > 0. ■ 

Example 13: Consider the process z(t) — X\e tt( - + X2e ltTI where f , rj £ R, £ ^ r\ and X\, Xi are 
proper, independent zero-mean Gaussians, i.e. ^X\Xi = EX1X2 = EX 2 = EXf = 0. The process 
z is proper, and WSS because the covariance function is r z (t,s) = E|Xi| 2 e^'* _s - ) + K\X2\ 2 e" 1 ^ t ~ s \ 
The real part is Re r z {t,s) = 2r x (t,s) = 2p x (t - s) = E\Xi\ 2 cos(£(t - s)) + E|X 2 | 2 cos^i - s)). 
From Proposition [T2] we may conclude that T = 0, that is, the IF process <p(t) has infinite variance 
for all t £ K, So a linear combination of pure exponential functions, with independent proper Gaussian 
weights, has an infinite-variance IF process everywhere unless it consists of a single term, in which case 
it follows from Remark [6] that its IF process is deterministic (with variance zero for all time points). 

Next we study the class of harmonizable processes which is larger than the class of mean-square 
continuous WSS processes. By the following three examples we show that T = R or T = may occur 
for harmonizable processes that are not WSS. We do not know whether f)CTCR may occur. 

Example 14: Let x(i) and y(t) be independent Gaussian real-valued processes with identical locally 
stationary Q, lfl8ll covariance function 

r x {t,s) = r y {t,s) = cxp(-2a(-^-) - ^ ~ ^ ) ' 

which is a covariance function provided b > a > [18|. We compute did2r x (t,t) = ((b — a) + 
4a 2 i 2 )exp(— 2at 2 ) and d\r x (t,t) = — 2aiexp(— 2at 2 ) which gives 

r x (t, i)<9i d 2 r x (t, t)-{d x r x {t, t)f = (b-a+4a 2 t 2 -4a 2 t 2 ) cxp(-4ai 2 ) = (b-a) exp(-4a£ 2 ) > 0, t £ 1 



if b > a. Thus by d52]l T = in this case. 

Example 15: Let g £ C 1 (R) be real-valued and the Fourier transform of a bounded measure, and 
let again x(t) and y(t) be independent Gaussian real-valued processes with equal covariance function 

r x (t, s) = r y (t, s) = g(t)g(s). If z{t) = x{t) + iy(t) then r z (t, s) = 2g(t)g{s). We have 

r x (t, t^dir^t, t) - (d x r x {t, t)) 2 = g(t) 2 g(t) 2 - (g(t)g(t)) 2 =0, t £ R, 

which means that T = R where T is defined by d52t . In this example \r z (t, s)\ 2 = r z (t,t)r z (s,s), i.e. 
we have equality in Cauchy-Schwarz for all t, s £ R, which implies that there exists so £ R such that 
z(t) — c(t)z(so) where c(t) is a deterministic function. In fact c(t) is a multiple of git). Thus we are 
in the situation of Remark [6] 

Example 16: This example is a generalization of Example [T3l Let z(t) — Xie 1 ^ + X2e ltri where 
£,77 £ R, £ 7^ 77 and Xl, X2 are jointly proper, zero-mean, unit variance Gaussians, i.e. E|Xl| 2 = 
E|X 2 | 2 = 1 and EX X X 2 = EX 2 = EX 2 = 0. In contrast to Example [Til we allow a non-zero 
correlation between X\ and X2, EX1X2 := c ^ 0. We assume \c\ < 1 since otherwise we have equality 
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IEX1X2I 2 = E|Xl| 2 E|X2| 2 in the Cauchy-Schwarz inequality, which implies that X2 is a complex 
multiple of X\, reducing the problem to that of Remark [6] 

The process z is proper, harmonizable but not WSS because the covariance function is r z (t,s) = 
e «(*-«) + ^n(t-s) + Ce m- V s) +ge »(i7t-M. The terms in the decomposition r z = 2r x + 2ir yx are 

2r x (t, s) = cos(£(i — s)) + cos(i](t ~ s)) 

+ Re c (cos(£i — rjs) + cos(r]t — £s)) + Im c (sin(??i — £s) — sin(£i — rjs)) , 

2r yx (t, s) = sin(£(7; - s)) + sin(r/(i - a)) 

+ Re c (sin(£r; — 77s) + sin(ryt — £s)) + Im c (cos(£t — 77s) — cos(r/t — ^s)) . 
Straightforward computations give 

2r x (t,t) = 2(1+ Re ccos(t(£- 77)) - Im csin(t(^- 77))), 

2<9i<Vx(M) = f 2 + ?7 2 + 2£tj (Re ccos(t(e - r?)) - Im csin(i(^ - 77))) , 

2d t r x (t,t) = (v-0 (Re csin(i(£ - 77)) +Im ccos(t(e - 77))) , 

2d iryx (t,t) = (77 + (1 + Re ccos(i(e - 77)) - Im csin(i(£ - 77))) , 
which in turn yields 

4 (r x (t, t)d x d 2 r x (t, t) - (d iryx (t, t)f - (d irx (t, t)) 2 ) = (£ - 77) 2 (1 - |c| 2 ) > 

for all tel. This means that T = 0. We conclude that a linear combination of two pure exponential 
functions, with correlated proper Gaussian weights of equal power, is a nonstationary process whose IF 
process ip(t) has infinite variance everywhere. 

VIII. Conclusion 

We have derived the pdf of the IF of a mean-square differentiable proper Gaussian stochastic process 
at a fixed and arbitrary time point. The proof does not involve the frequency domain. It shows that the 
IF has either infinite or zero variance. Moreover, if we assume that the process is harmonizable then 
we obtain as a consequence of our main result the formula for the mean of the IF as a normalized 
first-order frequency moment of the Wigner spectrum, again for fixed time. This is a generalization of 
the corresponding well-known formula for deterministic functions. 

The question of the dichotomy zero/infinite variance of the IF process is left somewhat open in this 
paper. We do not know if one process may exhibit nonempty time sets of both kinds. We have shown 
by examples that IF processes with constant zero variance exist, and IF processes with constant infinite 
variance exist, and each of these cases can occur both inside and outside the class of WSS processes. 
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